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Magnetodynamo Lifetimes for Rocky, Earth-Mass Exoplanets 
with Contrasting Mantle Convection Regimes 
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Abstract. We used a thermal model of an iron core to calculate magnetodynamo evo- 
lution in Earth-mass rocky planets to determine the sensitivity of dynamo lifetime and 
intensity to planets with different mantle tectonic regimes, surface temperatures, and core 
properties. The heat flow at the core-mantle boundary (CMB) is derived from numer- 
ical models of mantle convection with a viscous/pseudo-plastic rheology that captures 
the phenomenology of plate-like tectonics. Our thermal evolution models predict a long- 
lived (~8 Gyr) field for Earth and similar dynamo evolution for Earth-mass exoplanets 
with plate tectonics. Both elevated surface temperature and pressure-dependent man- 
tle viscosity reduce the CMB heat flow but produce only slightly longer-lived dynamos 
(~8-9.5 Gyr). Single-plate ("stagnant lid") planets with relatively low CMB heat flow 
produce long-lived (^10.5 Gyr) dynamos. These weaker dynamos can cease for several 
billions of years and subsequently reactivate due to the additional entropy production 
associated with inner core growth, a possible explanation for the absence of a magnetic 
field on present-day Venus. We also show that dynamo operation is sensitive to the ini- 
tial temperature, size, and solidus of a planet's core. These dependencies would severely 
challenge any attempt to distinguish exoplanets with plate tectonics and stagnant lids 
based on the presence or absence of a magnetic field. 



1. Introduction 

Many planets in the Solar System have or had magnetic 
fields produced by internal dynamos and this is presumably 
also true among planets orbiting other stars. Magnetic fields 
of exoplanets have yet to be unambiguously detected, but 
corotating regions of elevated chromospheric emission on the 
host stars of some close-in giant planets can be explained by 
magnetic interaction [Shkolmk et al., 2005; Shkolnik et al., 
2008; Fares et al., 2012] and one explanation for the inflated 
radii of many giant planets on close-in orbits is heating by 
ohmic dissipation during circulation of a partially ionized 
upper atmosphere through the planet's magnetic field [Baty- 
gin and Stevenson, 2010]. Past or present operation of mag- 
netodynamos has been linked to the interior structure and 
thermal evolution of the rocky planets of the inner Solar 
System [Breuer et al., 2010; Stevenson, 2010] and the detec- 
tion of magnetic fields might, in principle, provide insight 
into the interiors of exoplanets. 

In differentiated silicate-metal planets like the Earth, dy- 
namos are produced by convection within the liquid or par- 
tially liquid iron core if the magnetic Reynolds number 
Re m > 40 [Stevenson, 2010]. Dynamo operation is asso- 
ciated with ohmic dissipation, which requires entropy pro- 
duction in excess of that from conduction along the core's 
adiabat. Entropy production can be due to cooling of the 
core in excess of conduction, release of latent heat during 
solidification of an inner core, or formation of buoyant fluid 
generated by the exclusion of lighter elements such as sul- 
fur when iron crystallizes. All three sources require removal 
of heat from the core and thus core cooling is considered 
a requirement for dynamo operation. The decay of the 
long-lived radioactive isotope 40 K is another possible energy 
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source [Lee and Jeanloz, 2003; Murthy et al., 2003; Nimrao 
et al., 2004] although little such potassium may partition 
into a core [Corgne et al., 2007]. Planetary rotation rate 
is thought to have little influence on dynamo strength al- 
though it may control the relative strength of the dipole 
and multipolar components [Zuluaga and Cuartas, 2012]. 

In an Earth-like planet, the core is surrounded by a sili- 
cate mantle and heat must be carried from the core-mantle 
boundary (CMB) by subsolidus convection. The viscosity 
of a silicate mantle is strongly temperature-dependent and 
most likely the mantle contains long-lived radionuclides that 
provide heat and sustain high mantle temperatures over bil- 
lions of years. The ability of the mantle to carry heat from 
the core depends also on the style of mantle convection, dis- 
tinguished by the degree to which the brittle lithosphere 
participates in convection, e.g. plate tectonics or stagnant 
lid modes. Thus, magnetodynamo operation is tied to the 
tectonics and thermal evolution of the planet's mantle. The 
lack of a dynamo in Venus has been ascribed to the appar- 
ent absence of plate tectonics [Buffett et al., 1996; Nimmo, 
2002]. 

Previous theoretical studies have predicted the evolution 
of dynamos in rocky exoplanets [Gaidos et al, 2010; Driscoll 
and Olson, 2011; Tachinami et al., 2011]. Gaidos et al. 
[2010] examined the effect of surface temperatures and plate 
tectonics. They found that higher surface temperatures can 
enhance the dynamo due to the lower viscosity of a hot- 
ter mantle and that dynamos are weaker but still operate 
on planets with stagnant lids. To explain the absence of a 
present-day dynamo on Venus it was necessary to invoke an 
order-of-magnitude enhancement in mantle viscosity. They 
used parametric models to describe mantle convection and 
heat transport across the CMB. However, accounting for 
the complexity of subsolidus convection in a planetary man- 
tle with a strongly temperature-dependent viscosity requires 
numerical simulations. 

Furthermore, numerical models can adequately treat the 
effects of a non-uniform surface temperature. Rocky plan- 
ets very close to their parent star will be tidally locked, 
most likely with a 1:1 spin-orbit ratio, and may lack an at- 
mosphere due to erosion by the stellar wind [Khodachenko 
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et al, 2007]. In this case, temperatures will vary markedly 
between the irradiated and dark sides, and possible pat- 
terns of mantle convection include a degree-one (asymmet- 
ric) mode [Van Summeren et al., 2011]. This could enhance 
heat flow inhomogeneity across the CMB and suppress a 
dynamo if peak-to-peak variations are of the same order as 
mean values [Olson and Christensen, 2002]. 

Finally, the rheology of planetary mantles is expected to 
be pressure-dependent as well as temperature-dependent, al- 
though the pressure- dependence for super-Earths is a topic 
of controversy. If viscosity increases with pressure [Sta- 
menkovic et al., 2011], this will suppress convective vigor 
in the deep mantles of massive planets [Wagner et al, 2011; 
Tackley et al., 2013], which in turn will reduce heat flow 
across the CMB. Tachinami et al. [2011] predict that, above 
a critical planetary mass, enhanced viscosity terminates dy- 
namo operation. On the other hand, Karato [2011] argues 
for a decrease in viscosity with pressure by 2 to 3 orders of 
magnitude due to changes in the diffusion mechanism (from 
vacancy to interstitial) and a series of chemical and elec- 
tronic phase transitions over 0.1-1 TPa. If this is correct, 
the deep mantles of super-Earths will be less viscous and 
layered convection may occur. 

To more thoroughly investigate the conditions for the ap- 
pearance of dynamos in Earth-like exoplanets, we investi- 
gated the influence of differences in convective regime, sur- 
face temperature, pressure-dependence of mantle viscosity, 
and a range of core properties (radius, temperature, thermal 
conductivity, and solidus) on the lifetimes and intensities of 
magnetodynamos in planets with the same mass and interior 
structure as Earth. We show how dynamos can be longer- 
lived for planets with immobile surfaces compared to those 
with plate tectonics. We then determined the sensitivity of 
dynamo behavior to the planetary and physical parameters 
of our simulations, and assessed how informative dynamo- 
generated fields may be for probing the interior dynamics of 
exoplanets. 

2. Methods 

2.1. General modeling approach 

To explore parameter space, we calculated a series of 
dynamo evolution calculations for different mantle tectonic 
regimes, surface temperatures, and core properties. A cor- 
responding number of forward calculations for mantle evo- 
lution over billions of years is computationally impractical. 
Therefore, we adopted the following approach: (i) We cal- 
culated the steady state CMB heat flow, Qcmb, from nu- 
merical models of mantle convection with an imposed con- 
stant CMB temperature, Tcmb , and rate of mantle internal 
heating, H m . (ii) We repeated this procedure for different 
CMB temperatures to derive parameterizations of Qcmb 
vs. Tcmb for each convection scenario, (iii) We employed 
these parameterizations in an analytical model for the ther- 
mal evolution of an iron core, whi ch we s olved numerically. 
In these models, we rescaled the Qcmb values to account 
for time-dependent mantle internal heating rates, consistent 
with boundary layer theory, (iv) We calculated dynamo in- 
tensities using a scaling law with core entropy terms. 

2.2. Thermal model of the core 

We consider liquid/solid Fe cores with the same total ra- 
dius of that of Earth, except in Section 3.4 where we vary the 
core radius. To estimate the intensity of a magnetic dipole 
field at the CMB (Bcmb) we adopt the following scaling law 
derived from numerical dynamo simulations [Christensen 
and Aubert, 2006; Olson and Christensen, 2006; Aubert et 
al, 2009]: 



Bcmb — a\fjxop 



<I>T 

Wc (r C MB-n) 



1/3 



with $ the entropy production rate associated with ohmic 
dissipation due to electrically charged currents in a conduct- 
ing fluid, /xo the magnetic diffusivity, p the average core 
density, T the temperature at which ohmic dissipation oc- 
curs [Roberts et al., 2003], M c the core mass, and tcmb 
and ri the radii of the CMB and inner-outer core bound- 
ary, respectively. Parameters and values are listed in Table 
1. The dimensionless coefficient a — Ci ■ bdi P = 0.165 com- 
bines a prefactor used to extrapolate dynamo modeling re- 
sults from comparatively small Rayleigh numbers and large 
Ekman and magnetic Prandtl numbers to planetary condi- 
tions (ci = 1.65, Aubert et al. [2009]) and a fitting coefficient 
that relates the mean field strength inside the core to the 
dipole field at the CMB (b dlp « 10, Aubert et al. [2009]). 
The adopted scaling law successfully predicts the magnetic 
moments of the Solar System planets and Ganymede [Olson 
and Christensen, 2006]. For pure dipoles, the average field 
intensities at the planetary surface and CMB are related 
through B srfc « (r C MB/r sr f c ) 3 BcMB. 

We calculated the entropy production rate associated 
with ohmic dissipation, $, from the time-dependent entropy 
balance for the core: 



& + E K = Es + E L + E G . 



(2) 



The right-hand side of Eq. (2) contains the source terms 
associated with secular cooling (Es), release of latent heat 
(El), and release of gravitational potential energy (Eg). 
The sink term Ek is the entropy associated with maintain- 
ing thermal conduction along an adiabatic temperature gra- 
dient. A convection-driven dynamo requires $ > 0. Ex- 
pressions for the terms of Eq. (2) are listed in Table 2 and 
are based on Nimmo [2009] and Gaidos et al. [2010] and 
references therein. 

For a vigorously-convecting core one can assume a hy- 
drostatic, isentropic, and well-mixed state. The associated 
adiabatic density and temperature profiles can be approx- 
imated as (Labrosse et al. [2001]; Labrosse [2003]; Nimmo 
[2007]): 



and 



p(r) = pcmb exp 



T ad (r) — T C mb exp 



r 2 C MB - r 2 



L 2 



»(*)• 



D 2 



(3) 



(4) 



where pcmb and Tcmb are the CMB density and tempera- 
ture, respectively, and L and D are the density and thermal 
length scales, respectively (see Table 1). The assumed ther- 
mal equilibration between inner and outer core is reasonable 
if the time scale of inner core convection is shorter than that 
of inner core growth [Labrosse et al., 2001], a likely condition 
during the early stages of inner core growth [Buffett, 2009]. 
We calculated core solidification by expressing the solidus 
(T so i) as a Lindemann law of melting (e.g., Labrosse et al. 
[2001]) 



dlogp 



(5) 



where 7 is the Griineisen parameter. We ignored the possi- 
ble effect of unknown concentrations of light elements such 
as sulfur on the phase diagram of the core. In all our calcula- 
tions, the core solidifies from the planetary center outward. 
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Solidification at the top of the core (iron snow) is more likely 
for planets of two Earth masses or more [Gaidos et al, 2010]. 
We used the following energy balance equation for the 
core: 



Qcmb — Qs + Ql + Qg- 



(6) 



The source terms are the secular core cooling Qs, the la- 
tent heat release Ql, and the gravitational energy release 
Qg (specified in Table 2). The CMB heat flow Qcmb is the 
only relevant sink term and implicitly includes conductive 
heat transport in the core along an adiabat. We neglected 
core heating from the decay of radioactive isotopes such as 
40 K, as well as viscous dissipation [Braginsky and Roberts, 
1995]. Because Qs, Ql, and Qg all depend linearly on the 
cooling rate dTcMs/dt (Table 2), this cooling rate is readily 
solved if the heat flow Qcmb is known. We use dTcMB/dt 
to solve the entropy terms (Eq. 2, Table 2) and to update 
the core temperature profile (Eq. 4) and inner core radius 
at model time steps of 1 Myr. 

We characterized magnetodynamo evolution by (i) the 
onset time of inner solid core formation, t ons , when Ql and 
Qg become nonzero and the dynamo intensifies, (ii) the dy- 
namo cessation time, t mag , defined as the time at which the 
dipole field intensity reduces to zero due to core freeze-out, 
and (iii) the maximum dipole field intensity from the start of 
the model until core solidification is complete, B^^b- We 
neglected the possibility that part of the liquid core would 
not contribute to the magnetodynamo due to stable ther- 
mal and/or chemical stratification. However, we identified 
epochs of partial stabilization by calculating the onset times 
of thermal stratification, t 3tr at, i.e. when Qcmb < Qk with 
Qk the thermal conduction along the adiabatic temperature 
gradient at the top of the core (Table 2). For Earth, a strat- 
ified layer is probably <100 km thick [Gubbins, 2007] and 
its influence on the magnetodynamo is likely small. 

2.3. Dynamical model of the mantle 



We calculated Qcmb values from numerical models of 
mantle convection with constant internal heat production 
(4x 10~ 12 W kg -1 ). The methodology of our mantle convec- 
tion models has been described in detail in Van Summeren 
et al. [2011] and here we discuss only the aspects relevant 
for this paper. Mantle flow was described by a Boussinesq 
formulation of an incompressible fluid at infinite Prandtl 
number. The associated conservation equations for mass, 
momentum, and energy were solved using the finite element 
package CitComS-3.1.1 [Zhong et al, 2000]. The mantle do- 
main is modeled as a 360° annulus (Figure 1), with 513 x 65 
mesh nodes in the lateral and vertical direction respectively, 
and mesh refinement toward the top and bottom boundaries 
where the largest viscosity contrasts occur. 

To permit plate-like behavior of the surface boundary 
layer, we used a composite viscous/pseudo-plastic rheology 
model [Tackley, 2000]. In low-stress regions, we assigned an 
Arrhenius-type rheology: 



r}v(T') = r; exp 



23.03 



1 + T' 



(7) 



with reference viscosity r/o = 5 x 10 20 Pa s. The dimen- 
sionlcss potential temperature T' v scales to the dimensional 
potential temperature T p as T' v = (T p - 273)/2400. The 
viscosity r\ v decreases by 5 orders of magnitude when T v in- 
creases from 273 to 2673 K. In regions where the convective 
stress exceeds an assigned yield stress (a y ), we implemented 
pseudo-plastic yielding by recalculating of the viscosity r\ y = 
(j y /2tii, where eh is the second invariant of the strain rate 
tensor. This yielding allows for lithospheric break-up and 
concentrates strain in narrow regions. We computed the 



composite mantle viscosity as r\ m = min (r) v (T'),r) y ). We 
assigned a y = 150 MPa in our convection experiments, ex- 
cept for the stagnant lid case where pseudo-plastic yielding is 
precluded. This parameterization underestimates the tem- 
perature dependence of viscosity compared to experimen- 
tal flow laws for olivine (e.g. [Hirth and Kohlstedt, 2003]), 
but it avoids sharp viscosity gradients that are computation- 
ally expensive. This viscosity parameterization is commonly 
used in modeling the dynamics of planetary mantles and suc- 
cessfully reproduces the phenomenology of plate- like tecton- 
ics by localizing strain in regions that resemble subduction 
zones and spreading ridges [Tackley, 2000; Van Heck and 
Tackley, 2008; Foley and Becker, 2009; Van Summeren et 
al, 2011]. Although the viscous strength of plates is under- 
estimated, plate tectonics is successfully reproduced because 
surface deformation occurs by yielding rather than viscous 
deformation. 

For each scenario, we computed 8 cases with CMB po- 
tential temperatures in the range Tcmb p = 1873-3473 K 
(Tomb,, =0.67-1.33). We calculated physical (or "real") 
CMB temperatures (Tcmb) from potential temperatures 
(Tcmb p ) by correcting for the influence of adiabatic com- 
pressibility: i.e., using Tcmb = Tcmb p exp(Di), with dis- 
sipation number Di — a m gh/cp m — 0.46, mantle thermal 
expansivity a m = 2 x 10~° K _1 , mantle gravitational ac- 
celeration g — 10 m s~ 2 , mantle thickness h — 2.891 x 10 6 
m, and mantle heat expansivity cp m = 1250 J kg -1 K _1 . 
Because each model has a constant CMB temperature and 
heat production, Qcmb e ventually reaches a statistically 
steady state value, Qcmb- Typically, we ran each case for 
a model time of ~9 Gyr (i.e. 32,000 model time steps) and 
we removed the first ~5 Gyr (12,000 time steps) to exclude 
initial transients. For the results presented in Section 3.4 we 
increased the maximum core temperature to Tcmb — 5000 
K by extrapolating the heat flow data, using linear regres- 
sion to the filled symbols shown in Figure 3b. We applied 
the same technique to estimate Qcmb for Tcmb >4500 K 
for scenario VI (i.e. the case represented by the right-hand 
brown circle in Figure 3a relies on extrapolation). 

To account for a time-dependent rate of mantle internal 
heat production, we rescaled the CMB heat flow using the 
quasi steady state condition QcMB{i) = Q S rfc{t) — H m (t). 
We scaled the surface heat flow by Q sr fc(t) oc Ra H (i), with 
Rau{t) = a m PmgH m (t)h 3 /k m K m rj m the Rayleigh number 
for an internally heated fluid, fc m =5.65 W _1 K the man- 
tle thermal conductivity, and p m = 4500 kg m~ 3 the man- 
tle reference density. We thus have Q ar fc(t) oc H m {t). 
We also assumed a quasi steady state heat balance for the 
unsealed terms (indicated with subscripts "0"): Q S rfc — 
Qcmb + Hm ■ We then rescaled the CMB heat flow with 



j(i) = (Qcmb + H mo ) 



H m (t) 



1/3 



tin 



ti m (t) (8) 



We used /3 = 0.25, consistent with boundary layer theory 
for an internally heated mantle. Deviations from the theory 
may arise, for example due to contributions from bottom 
heating or the dependence of Ran on mantle heating rate 
through rim- We validated that /3=0.25 is appropriate for 
our models by calculating the surface heat flow from two 
verification experiments with mantle internal heating rates 
2 x 10~ 12 and 6 x 10~ 12 W kg -1 (Figure 2), which yielded 
j3 = 0.27. We evolved the mantle heat production rate with 
respect to time as: 



H m (t) = Y, ClWex P 



-tin 2 



'1/2 



(9) 



with isotope concentration C, heat production rate H, and 
half-life time t 1/2 , for radiogenic isotopes 238 U, 235 U, 232 Th, 
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and 40 K (values from Turcotte and Schubert [2002], Ta- 
ble 4-2). Estimates for present-day mantle heating (e.g. 
7.4 x 10~ 12 W kg~\ [Turcotte and Schubert, 2002]) depend 
on assumptions of the Urey number. To account for this 
uncertainty, we considered a range of H mo values (lxl0~ 12 
- 8xl0 -12 W kg -1 ) and adjusted the heat production rates 
for radiogenic isotopes H l (Eq. 9), accordingly. 

For comparison with boundary layer theory, we recast 
our results in terms of dimensionless parameters and express 
Qcmb as 

QCMB — 4:TTrcMBk m (TcMB p — < T mp >)NUQMB, (10) 

with Nucmb the Nusselt number for the bottom boundary 
layer. Boundary layer theory predicts Nucmb oc Ra,fl MB 
for a bottom-heated mantle [Schubert et al., 2001]. We de- 
fine the Rayleigh number for the bottom boundary layer as 



RacA 



p m ga m (TcMB p — < T mv >)h 3 



(11) 



2.4. Parameter values and model scenarios 



We calculated CMB heat flows for 5 scenarios of mantle 
convection and temperature (Table 3). We simulated a nom- 
inal Earth ( "cold" or C model) with a uniform surface tem- 
perature T sr fc — 273 K. Using a "hot" (H) model with a uni- 
form T sr f c — 759 K we simulated a close-in planet (a=0.13 
AU) with efficient surface heat redistribution. An "asym- 
metric" (A) model was used for close-in planets (a=0.13 AU) 
that lack heat redistribution, and we assigned a surface tem- 
perature that is elevated at the substellar point (T sr / C =1073 
K), decreases gradually to the terminus, and is constant on 
the night side (T sr f c — 273 K). To investigate planets with 
an immobile surface and convecting interior, such as Venus, 
we ran a "stagnant lid" (SL) model by removing pseudo- 
plastic behavior from the viscosity parameterization. For 
the "viscosity increase" (VI) model we assigned an expo- 
nential increase of the viscosity by a factor of 10 across the 
mantle depth range, in addition to temperature-dependent 
viscosity (Eq. 7). 

We complemented our investigation of different convec- 
tion regimes (Sections 3.1-3.3) with a study of the sensi- 
tivity of results to core properties (Section 3.4). We cal- 
culated the sensitivity of the dynamo cessation time (J m „ s ) 
and maximum CMB field intensity {Bq^b) 1° t wo plane- 
tary properties: core radius (rcMB) and initial CMB tem- 
peratures (T^mb ) , all for the nominal C scenario of mantle 
tectonics We also calculated the sensitivity of our results to 
core solidus (T so i) and heat conductivity (fc c ), parameters for 
which the values at pressures corresponding to Earth's cen- 
ter are controversial [Morard et al, 2011]. Our solidus range 
equates to a range of the Griineisen parameter 7 € [1-1- 
1.4], reflecting current uncertainties in high- pressure mineral 
physics. Where we varied rcMB, we recalculated the inte- 
rior structure using the interior model described in Gaidos 
et al. [2010]. This uses the third-order Birch-Murnagham 
(BM) equations of state, and includes the Thomas-Fermi- 
Dirac contribution to the pressure using the formulation of 
Zapolsky and Salpeter [1969]. The density in the liquid part 
of the core is adjusted by a fixed fraction Sp/p to account 
for the presence of light elements [Lee and Jeanloz, 2003]. 
We adjusted the solidus at the planetary center according 
to Eq. (5). For a valid comparison we adjusted the mantle 
thickness to maintain a constant Rayleigh number, i.e. the 
product of mantle density, gravitational acceleration, and 
the cube of the mantle thickness, evaluated at the CMB. 

3. Results 

3.1. Mantle dynamics and CMB heat flow 



Figure 3a shows our computed CMB heat flow Qcmb as 
a function of Tcmb for the 5 mantle convection scenarios 



described in Section 2.3 and listed in Table 3. Models C, 
H, and A build on the work of Van Summeren et al. [2011] 
and their dynamics will be discussed only briefly. The nom- 
inal Earth-like scenario C with a cold surface (T' sr / C =273 
K) is characterized by plate-like behavior with rigid litho- 
spheric plates that are broken by narrow regions of high 
strain that resemble slab-like downwellings and ridge-like 
upwellings (Figure 1). The Qcmb values increase with in- 
creasing Tcmb (Figure 3a, blue curve), which reflects ele- 
vated heat transport in a more vigorously convecting mantle. 

In scenario H of a close-in planet with a uniformly 
elevated surface temperature (T sr / C =759 K), the surface 
boundary is too weak to maintain coherent plates. Instead, 
the surface material is highly mobile and deforms diffusely, 
a possibility for close-in planets with surface temperatures 
close to the silicate solidus. We limited the temperature 
dependence of viscosity for numerical reasons and therefore 
likely underestimated the viscosity of the surface boundary 
layer. Planets with T sr f c — 759 K may instead have immobile 
or episodically overturning surfaces [Armann and Tackley, 
2012], and we investigated this stagnant lid behavior with 
our SL models, described below. For scenario H, the Qcmb 
values are only slightly smaller than for scenario C (Figure 
3a, cf. red and blue curves). Indeed, the elevated surface 
temperature causes a mantle temperature increase that re- 
duces the CMB heat flow, but this is mostly offset by a lower 
viscosity in the hot boundary layer that enhances CMB heat 
flow. 

In scenario A, the hemispheric contrast in surface temper- 
atures causes a tectonic dichotomy with diffuse deformation 
on the hot day side (similar to scenario H) and plate-like 
tectonics on the cold night side (similar to scenario C) [ Van 
Summeren et al., 2011]. This mixed style of surface defor- 
mation results in Qcmb values that are intermediate those 
of scenarios C and H (Figure 3a) . An asymmetric (harmonic 
degree 1) convection pattern develops in the mantle interior. 
This pattern involves convective upwellings that are concen- 
trated near the hot sub-stellar point, near-surface flow from 
the hot day side to the cold night side where downwellings 
concentrate, and a deep mantle return flow toward the day 
side. A small hemispheric heterogeneity in CMB heat flow 
of ~10% in scenario A is unlikely to destabilize a dynamo 
[Olson and Christensen, 2002]. 

In scenario VI (in which viscosity increases with pressure) 
CMB heat flows are ~15-25% lower than in scenario C (Fig- 
ure 3a). This reflects the depth- increasing viscosity which 
suppresses deep mantle convective overturn and hinders de- 
scending cold slabs from covering and cooling the core. As 
a result, a thicker, less conducting CMB boundary layer de- 
velops, relative to reference scenario C. 

The stagnant lid convection scenario SL has an immobile 
surface through which heat is conducted, but vigorous con- 
vection still occurs in the underlying mantle. The Qcmb 
values are lower than for models C, H, and A for the Tcmb 
range we investigated (Figure 3b). This is because slow heat 
transport through the immobile surface reduces the temper- 
ature difference between the mantle interior and the core. 

We expressed our results in dimensionless quantities, i.e. 
the Nusselt and Rayleigh numbers of the bottom boundary 
layer as defined by Eq. (10) and Eq. (11). When bot- 
tom heating dominates (i.e. for Urey number U <0.5), our 
results approach the theoretical 1/3 power-law relationship 
for a bottom-heated mantle [Schubert et al., 2001] (Figure 
3b). Calculated Nucmb values deviate from the theoretical 
value when internal heating dominates (U >0.5), which oc- 
curs for Tcmb <4000-4500 K in our models (Figure 3b). For 
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U >0.5, Nucmb is larger for scenarios H and A than for sce- 
nario C (Figure 3b) because (locally) elevated surface tem- 
peratures thermally buffer the mantle interior and reduce 
the temperature contrast with the core (Eq. 10). The CMB 
heat flow also decreases, but this is mostly compensated by 
more efficient heat transport resulting from decreased man- 
tle viscosity. Compared to scenario C, scenarios VI and SL 
have a lower Nucmb which reflects a lower CMB heat flow 
due to the weaker influence of cold slabs in the deep mantle. 

3.2. Nominal thermal evolution of the core 

We calibrated the nominal Earth-like scenario by repro- 
ducing Earth's present-day inner core radius of 1220 km and 
heat flow of 5-15 TW [Lay et al., 2008]) at a model time of 
4.5 Gyr. In calculating the thermal evolution of the core, 
we accounted for time-dependent mantle heating by rescal- 
ing the parameterized CMB heat flow values (Figure 3, blue 
curve) and investigated the range of Tcmb and H mg val- 
ues shown in Figure 4. Different combinations of parameter 
values can reproduce Earth's inner core radius. For exam- 
ple, the calculated present-day inner core radius increases 
for progressively lower Tcmb, lower H mo , or higher solidus 
T so i (i.e., larger Griineisen parameter 7) and the present- 
day Qcmb increases for higher Tcmb, lower H m „, or higher 
T so i. We chose T£$ s =5000 K, # mo =3.8x 1(T 12 W kg" 1 (la- 
bel 1 in Figure 4), and 7=1.3 (T s c 7=5391 K) and the corre- 
sponding core evolution (Figures 5a-c) serves as a reference 
for other scenarios. Dynamo operation starts at t ~0.8 Gyr 
(<f> > 0, Figure 5c), when the CMB heat flow has become 
sufficiently large (~14 TW) after an initial period of strong 
mantle heating that suppresses the CMB heat flow (Figure 
5b). After inner core formation (~4.1 Gyr), latent heat and 
gravitational energy release contribute to the entropy pro- 
duction in the core, which causes an abrupt increase in <f> 
(Figure 5b) and a stronger dynamo. As the core contin- 
ues to solidify, <f> gradually decreases due to the progressive 
thinning of the fluid iron layer. Finally, dynamo operation 
ceases at ~8.2 Gyr when core solidification is complete. The 
onset of inner core formation coincides with a change in the 
Qcmb slope from negative to nearly fiat (Figure 5b, ~4.1 
Gyr). This is due to (i) the generation of latent heat and 
gravitational potential energy, which partly replenishes heat 
carried from the core and maintains high core temperatures, 
and (ii) cooling of the mantle due to decreasing internal heat 
production, which acts to enhance Qcmb- 

We calculated the sensitivity of dynamo evolution to H mo 
and Tcmb using 3 combinations that reproduce Earth's 
present-day inner core radius (values as specified in Fig- 
ure 4). A progressive decrease of T^mb by 500 K lowers 
the CMB heat flow, which is reflected in a corresponding 
decrease in $ and a delay in the onset time of dynamo op- 
eration (Figures 4 and 5d, labels 1-3). Among the models 
that reproduce Earth's present-day inner core radius, only 
those with Tqj^ b > 5000 K correspond with an early non- 
magnetic epoch <1 Gyr (label 1 in Figures 4 and 5d), in 
agreement with evidence for Earth's long-lived (>3.5 Ga) 
magnetic field [Biggin et al., 2008]. We therefore prefer case 
1 and use the corresponding values in subsequent models, 
unless otherwise mentioned. 

3.3. Dynamo evolution for different mantle convection 
regimes 

Figure 6 shows the evolution of the dynamo in the 5 man- 
tle convection scenarios (Section 2.4, Table 3). Throughout 
most of the evolution, the CMB heat flows for scenarios with 
plate tectonics or diffuse surface deformation (C, H, A, VI) 
are larger by ~2-4 TW than for t he stagn ant lid scenario SL 
(Figure 6a), a consequence of the Qcmb vs. Tcmb system- 
atics (Figure 3) . The impact of these contrasting heat flows 
on dynamo evolution is shown in Figure 6b. For scenarios 



with active surface recycling (C, H, A, and SL), dynamo 
operation starts at ~0.8-1.6 Gyr and continues until core 
freeze-out (~8.2-9.5 Gyr). For the stagnant lid scenario SL, 
an early dynamo ceases at ~4.1 Gyr because of the low CMB 
heat flow and corresponding low entropy production in the 
core. After a magnetic quiet period (~4. 1-5.0 Gyr), the 
dynamo can restart due to inner core formation and contin- 
ues until ~10.6 Gyr, i.e. longer than for all other scenarios 
that we investigated. The longevity of this second-stage dy- 
namo is due to late freeze-out in a slowly cooling single-plate 
planet with relatively low CMB heat flow. Unlike the other 
scenarios, thermal stratification can occur during a substan- 
tial part of the evolution of this core (£ s trat=4.4 Gyr, Table 
3). The CMB dipole field strengths Bcmb are of similar 
magnitude in the 5 different scenarios due to the weak de- 
pendence of Bcmb on $ (Eq. 1) and this is reflected in the 
limited variation in the maximum field intensities in these 
scenarios (Bg^ =260-310 ^T, Table 3). 

3.4. Sensitivity of dynamo evolution to core properties 

We investigated the sensitivity of dynamo evolution to 
core radius tcmb, initial temperature T c n ^ B , and material 
properties (solidus T so i and thermal conductivity k c ) in the 
nominal Earth-like scenario. Figures 7a and b show that 
hotter and bigger cores prolong dynamo operation. Planets 
with Tcmb <3500 K are non-magnetic because they (unre- 
alistically) start with and retain a solid core. A progressive 
increase in Tcmb causes a subsequently smaller increase of 
dynamo cessation times (Figure 7a) because hotter cores 
lose their heat more efficiently (Figure 3a). 

Recent first principle computations of liquid iron mix- 
tures [De Koker et al, 2012; Pozzo et al, 2012] suggest a 
core thermal core conductivity 2-3 times higher than previ- 
ous estimates that relied on extrapolations, e.g. Stacey and 
Anderson [2001]. Our results show that dynamo cessation 
times t mag are insensitive to k c in the range 60-240 W m _1 
K _1 (Figures 7c and d). An increase in k c increases the en- 
tropy term related to heat conduction Ek and this weakens 
the dynamo (lowers $, Eq. 2). For k c >160 W m _1 K _1 , 
the dynamo collapses prior to inner core formation, although 
this does not affect the cessation times t mag (Section 2.2). 
Dynamo cessation after inner core formation requires even 
higher k c (>500 W m _1 K -1 , outside the range shown in 
Figure 7c) to offset the latent heat and gravitational energy 
release contributions. Dynamo operation is insensitive to fc c 
after inner core formation and for k c <500 W m _1 K^ 1 be- 
cause: (i) the time of core freeze-out is controlled by Qcmb , 
which depends on the convective vigor of the mantle but is 
independent of k c , (ii) the entropy sink term Ek reduces 
to zero and this removes the dependence of $ on k c (Ta- 
ble 2). These k c limits are derived using a central solidus 
T so f=5391 K- Dynamo lifetimes decrease for progressively 
larger T so i (Figures 7c and d) because planetary cores so- 
lidify at higher temperatures, i.e. generally earlier in the 
planet's history. 

Peak field intensities are only weakly sensitive to core 
properties (tcmb, Tcmb, T so i, k c ) because maximum in- 
tensities occur at inner core formation (Figures 5c and 6c) 
which occurs at a specific temperature regardless of the on- 
set time and because Bcmb ~ "I? 1 ' 3 (Eq. 1). Intensities 
decrease sharply, however, when dynamo lifetimes reduce to 
zero, i.e. for low initial core temperatures (T c r MB ^3500 
K) and high solidus (T s c e , n >8400, Figures 7b and d, solidus 
limit out of the range shown). 



4. Discussion 

4.1. Comparison with Earth and Venus 
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Our nominal Earth-like scenario successfully reproduces 
evidence for a geodynamo that has operated since at least 3.5 
Ga [Biggin et al., 2008], and Earth's present-day inner core 
radius of 1220 km. The corresponding present-day CMB 
temperature of 3875 K in our models is compatible with ex- 
perimental and computational estimates for iron alloys of 
3700-4300 K (e.g., Boehler et al. [2000]) and estimates of 
3600-4200 K based on seismologically inferred deep man- 
tle structures combined with calculated elastic properties of 
lower mantle minerals [Lay et at, 2006; Van der Hilst et al., 
2007; Kawai and Tsuchiya, 2009]. These constraints restrict 
the evolution of core energetics, which critically depends on 
the CMB heat flow. In our preferred scenario, the CMB 
heat flow is ~15 TW at £=4.5 Gyr, consistent with recent 
estimates of the present-day CMB heat flow (5-15 TW, Lay 
et al. [2008]). The CMB heat flow varies between ~14 and 
17 TW, except for the first ~1 Gyr when strong mantle 
heat production impedes core cooling. Our calculations do 
not require a source of radiogenic heat in the core to match 
Earth's long-lived dynamo, 

A relatively young inner core (~500 Ma) and a maximum 
rate of entropy production associated with ohmic dissipation 
before and after inner core formation of respectively $ ~100 
MW KT 1 and ~1000 MW K" 1 are consistent with other 
core evolution studies with comparable CMB heat flow val- 
ues [Labrosse, 2003; Nimmo, 2007]. Previous studies have 
computed older cores (~1 Ga), likely because of a lower 
Qcmb ~10 TW [Labrosse, 2003; Breuer et al., 2010; Gaidos 
et al., 2010]. We prefer relatively high values that prevent 
the formation of a stably stratified layer in the liquid core 
(requiring Qcmb > Qk ~ 11-14 TW, Figure 5b) which, for 
Earth, is absent or thin (<100 km, [Gubbms, 2007]). 

In this light, additional mechanisms that can modulate 
Earth's CMB heat flow are important. Enhanced CMB heat 
flow is possible from the decay of 40 K in Earth's core, and 
destabilization of the mantle bottom boundary layer due 
to a viscosity reduction within deep mantle post-perovskite 
regions [Nakagawa and Tackley, 2011]. It is possible that 
Earth's viscosity increase with depth is larger (factor of 
~1000, Mitrovica and Forte [2004]; Steinberger and Calder- 
wood [2005]) than employed in our models. Our results show 
that an increase of the pressure dependence reduces convec- 
tive overturn in the deep mantle, which lowers the CMB heat 
flow and weakens the dynamo (cf. models C and VI, Figure 
6), consistent with parameterized convection models [Tachi- 
nami et al., 2011]. Dynamo operation may be viable for 
planets with strong pressure-dependence, however, because 
of sustained convective heat loss due to feedback between 
viscosity, temperature, and internal heating [Tozer, 1965; 
Tackley et al., 2013]. Possible mechanisms that can decrease 
the CMB heat flow include a reduced radiative conductivity 
due to high-pressure transitions from high-spin to low-spin 
in iron atoms [Badro et al., 2004; Goncharov et al., 2006], 
thermal insulation due to a compositionally-dense layer in 
the deep mantle [Nakagawa and Tackley, 2004, 2010], and 
reduced convective heat transport due to a decrease of the 
thermal expansivity with pressure [Chopelas and Boehler, 
1992]. 

It has been widely debated why a dynamo is absent on 
Venus, which has a similar size and interior structure as 
Earth but lacks plate tectonics and has a higher surface 
temperature. Our stagnant lid model results raise the pos- 
sibility that a dynamo is currently inactive in Venus because 
of a low CMB heat flow in the absence of a crystallizing inner 
core. If true, a dynamo may restart after the onset of core 
freezing which enhances the entropy production. This sce- 
nario is consistent with previous calculations of Stevenson 
et al. [1983] and would not require additional mechanisms, 
such as reduced mantle convection due to dehydration stiff- 
ening [Nimmo and McKenzie, 1996; Nimmo, 2002; Gaidos 
et al., 2010]. Indeed, the influence of an intrinsic viscos- 
ity increase may be limited by strong viscosity-temperature 
feedback [Tozer, 1965]. 



4.2. Uncertainties and model approximations 

Several uncertainties and model approximations temper 
strong statements about dynamo operation in rocky planets. 
Our viscosity parameterization underestimates the tempera- 
ture dependence of viscosity of mantle silicates compared to 
laboratory experiments. Planetary mantles with a stronger 
temperature dependence likely experience a smaller temper- 
ature decrease with time because thermal convection tends 
to be more self-regulating [Tozer, 1965]. If this causes a 
smaller temperature difference with the core, a lower CMB 
heat flow is expected to weaken dynamo operation. On the 
other hand, the bottom thermal boundary layer would have 
a lower viscosity, which locally enhances the mobility of up- 
wellings and this could increase the CMB heat flow and 
produce a stronger, but shorter-lived dynamo. For planets 
with immobile surfaces, Li and Kiefer [2007] showed that a 
stronger temperature dependence results in a thicker upper 
boundary layer and this reduces the surface heat flow and 
increases mantle temperatures, which reduces the CMB heat 
flow. A reduced CMB heat flow likely delays inner core for- 
mation and hampers dynamo operation until an inner core 
forms and this strengthens our explanation for a present-day 
non-magnetic Venus. 

In rescaling the CMB heat flow to account for time- 
dependent mantle heating we have neglected secular cooling 
of the mantle. Secular mantle cooling would contribute to 
the energy balance in a similar way as the radiogenic mantle 
heat production. We may therefore overestimate the CMB 
heat flow and overestimate dynamo intensities although, for 
Earth, secular cooling is likely ~3 times smaller than in- 
ternal heating [Turcotte and Schubert, 2002]. However, a 
larger early CMB heat flow would develop if the early core 
was super-heated relative to the mantle during an out-of- 
cquilibrium thermal state immediately following core- mantle 
differentiation [Stevenson, 1990]. Also, magmatism, while 
not included in our models, can enhance the CMB heat 
flow by cooling the mantle and this may be particularly 
important during Earth's early hot stages [Nakagawa and 
Tackley, 2012] and would promote early dynamo operation. 
To avoid the formation of an unrealistically large inner core 
may require a super-heated core or other compensatory heat 
source. 

Several processes are not considered in our models for 
Earth-sized planets but may drive or modulate dynamos of 
Solar System bodies of different size and structure. Mantles 
much thinner than Earth's allow for efficient core cooling 
and the growth of an inner core that promotes a chemically- 
driven dynamo, as has been proposed for present Mercury 
[Stevenson et al, 1983]. In contrast, a relatively thick man- 
tle and small core may hamper dynamo operation; this has 
been proposed to explain the absence of a long-lived dy- 
namo on the Moon [Runcorn et al., 1975]. Although Mars 
is less massive than Venus, early dynamo cessation due to 
insufficient CMB heat flow is also a possibility, perhaps trig- 
gered by a transition from active tectonics to stagnant lid 
convection [Nimmo and Stevenson, 2000]. Giant impacts 
can amplify or generate magnetic fields, which may have 
contributed to the Moon's paleomagnetic field [Hood and 
Vickery, 1984]. Mechanical stirring arising from mantle-core 
differential motion is another possibility for dynamo oper- 
ation on the early Moon and large asteroids [Dwyer et al., 
2011]. Stable stratification in the outer region of a liquid 
core, with a dynamo operating only at depth, can buffer a 
magnetic field and this may cause the weak field of Mercury 
[Christensen, 2006]. A chemically driven dynamo powered 
through rise or fall of (sulfur-rich) iron snow has been sug- 
gested for Ganymede [Bland et al., 2008; Hauck et al., 2006]. 
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First-principle computations of liquid iron mixtures [De 
Koker et al., 2012; Pozzo et al., 2012] suggest a thermal con- 
ductivity that is 2-3 times higher than previously derived 
values from extrapolations (k c ~60 W _1 m -1 , e.g. Stacey 
and Anderson [2001]). We performed the nominal Earth- 
like scenario at k c =120 W~ m _1 K _1 . We demonstrated 
that in the 60-240 W m _1 K _1 range, k c has no influence 
on the timing of dynamo cessation due to core freeze-out. 
However, k c >160 W m _1 K _1 may prevent dynamo oper- 
ation before inner core formation, which would contradict 
evidence for a terrestrial magnetic field since at least ~3.5 
Ga. For our stagnant lid scenario SL, lower k c would likely 
shorten the magnetic quiet period. This highlights the im- 
portance of at least establishing upper limits to the thermal 
conductivity under core conditions. 

4.3. Implications for exoplanets 

Our results suggest that Earth-mass rocky exoplanets 
with active surface deformation (plate tectonics or diffuse 
surface deformation) can have qualitatively similar dynamo 
evolutions and operate for ~8.2-9.5 billions of years, irre- 
spective of surface temperature (273-759 K). A longer- lived 
(~10.5 Gyr) dynamo is possible with stagnant lid convec- 
tion, due to the lower CMB heat flow. These dynamo life- 
times raise the question whether the interior dynamics of 
Earth-mass rocky exoplanets can be inferred by combining 
magnetic field detections and planet age estimates. The ages 
of stars (and their planets) can be inferred with a precision 
that would be sufficient for identifying planets in different 
phases of core thermal evolution, i.e. ~l-2 Gyr [Mama- 
jek and Hillenbrand, 2008]. Unfortunately, our results sug- 
gest that uncertainties in planetary rock properties make 
a distinction between planets with active surface recycling 
and stagnant lid planets extremely challenging. We demon- 
strated that characterization of exoplanet interior dynamics 
requires well-determined radius, initial temperature, solidus, 
and thermal conductivity of the core (Figure 7). In particu- 
lar, the strong sensitivity to the solidus emphasizes the need 
for constraining the equation of state of iron under core pres- 
sures and temperatures. 

What dynamo behavior can be expected for planets more 
massive than Earth? For rocky planets larger than ~2 Earth 
masses, it takes >4.5 Gyr to form an inner core and dynamo 
action before that time must be maintained by thermal con- 
vection [Gaidos et al, 2010; Driscoll and Olson, 2011]. The 
pressure-dependence of mantle viscosity is still debated [Sta- 
menkovic et al, 2011; Karato, 2011] but has important im- 
plications for mantle dynamics. Our results demonstrate 
that a viscosity increase of one order of magnitude across 
the mantle depth will hamper deep mantle convection and 
attenuate but not stop dynamo operation. Because terres- 
trial planets more massive than Earth have higher mantle 
pressures, their material properties probably depart more 
radically from terrestrial values. This makes it even more 
challenging to reliably predict dynamo lifetimes for massive 
Super-Earths. 

4.4. Detection of magnetic fields on exoplanets? 

Charged particle belts in the magnetospheres of Earth 
and the giant planets emit electron cyclotron emission at 
MHz frequencies. Magnetodynamo-hosting planets around 
other stars could presumably emit at similar frequencies and 
power levels [Hess and Zarka, 2011]. Searches for Jupiter- 
like emission from giant exoplanets have been undertaken 
[Bastian et al, 2000; Lazio and Farrell, 2007; George and 
Stevens, 2007; Farrell et al, 2004; Lazio et al, 2010a, b; 
Lecavelier Des Etangs et al, 2011]. However, predicted fre- 
quencies for Earth-like dynamos are below the 10 MHz cutoff 
of Earth's ionosphere. Moreover, an empirical relation be- 
tween emitted power and magnetic moment [Zarka, 2007] 
predicts that power emission is well below the detection 



threshold of even the most sensitive radio telescope (LO- 
FAR) [Farrell et al, 2004; Driscoll and Olson, 2011; Lazio 
et al., 2010b]. Detection would require a fortuitous combi- 
nation of increased stellar activity [Griefimeier et al, 2005], 
beamed emission in the direction of Earth, or intensification 
by a stellar flare [Driscoll and Olson, 2011]. Detection might 
also be possible with radio telescopes in space, which would 
not suffer from the low-frequency ionospheric cutoff. 

The magnetic field of a planet on a close-in orbit can in- 
teract with the field of its host star; periodic chromospheric 
emission from the host stars of some "hot Jupiters" could be 
the result of connectivity between the fields of the planet and 
the star [Shkolnik et al., 2008; Lanza, 2011]. Other observ- 
able phenomenae are associated with a planetary dynamo: 
A magnetic field could alter the circulation and tempera- 
ture distribution of a partially ionized upper atmosphere of 
a "hot" Earth [Castan and Menou, 2011], perhaps changing 
the pattern of reflected or emitted light, although these ef- 
fects have yet to be explored. A magnetized "tail" of gas es- 
caping from a planet could polarize transmitted light [Tachi- 
nami et al., 2011] or a bowshock could produce an asymme- 
try or variability in the in the primary occultation signal of 
a transiting planet [Vidotto et al., 2011a, b] 

The persistence of an atmosphere on a low-mass planet 
close to its host star is possible indirect evidence for a plan- 
etary dynamo, because a strong magnetic field protects the 
planet's atmosphere against erosion by a stellar wind and 
coronal mass ejections. The dense plasma environment close 
to the host star can remove hundreds of bars equivalent at- 
mosphere from an unmagnetized Earth-size planet over sev- 
eral Gyr [Khodachenko et al., 2007; Lamrner et al., 2007]. 
Any primary atmosphere dominated by H2 and He would 
be rapidly removed unless protected by a magnetic field. 
The presence of an atmosphere on a transiting planet can 
be detected by absorption features in spectra obtained dur- 
ing a transit [Charbonneau et al., 2002; Bean et al., 2010] or 
by its re-distribution of heat [Gaidos and Williams, 2004; 
Knutson et al., 2007]. Limitations to using planetary at- 
mospheres to infer a magnetodynamo are that planets can 
conceivably accrete without atmospheres [Raymond et al., 
2007; Lissauer, 2007] or may have lost their atmospheres 
early in their evolution due to increased stellar activity and 
before any dynamo started. We calculated an initial non- 
magnetic period of ~ 1—1.5 Gyr for all our model scenarios 
(Table 3) due to high rates of radiogenic mantle heating, 
which leaves all such planets vulnerable to loss of an atmo- 
sphere. Secondary atmospheres of volcanic CO2, H2O, and 
N2 may form but could also be removed by thermal escape 
and stellar wind erosion [Khodachenko et al., 2007; Tian, 
2009]. These possible complications will challenge attempts 
to infer the operation of a dynamo through the presence or 
absence of atmospheres on exoplanets. 
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Table 1. Planetary core parameters and values used in this 
study. Subscript conventions are "0" for zero pressure, "srfc" 
for surface, "m" for mantle, "c" for core, "cen" for planetary 
center, "CMB" for core-mantle boundary, and "i" for inner 
core. Values annotated with an asterisk (*) are specific for 
the nominal-like Earth model C (Table 3). Values are from 
(a) [Wang et al., 2002], (b) [De Koker et al., 2012; Pozzo et 
al., 2012], (c) [Anderson and Ahrens, 1994], (d) [Dziewon- 
ski and Anderson, 1981], (e) [Poirier and Shankland, 1993], 
(f) [Dziewonski and Anderson, 1981; Labrosse, 2003], and (g) 
[Buffett et al, 1996]. 
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Table 2. Terms of the energy (Q) and entropy (E) bal- 
ance (Section 2.2). Subscripts S, L, G, and K refer to secu- 
lar cooling, latent heat release, gravitational energy release, 
and conductive cooling, respectively. Other symbols are given 
in Table 1. The quantity F in the gravitational terms is 
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Table 3. Specification for the 5 mantle convection models 
C, H, A, SL, and VI (Section 2.4). Surface temperatures T sr f c 
are uniform across the planet, except for model A, where the 
surface temperature is 1073 K at the substellar point and de- 
creases sinusoidally toward the terminus and is kept constant 
on the night-side at 273 K. For model SL, pseudo-plastic rhe- 
ology is excluded [a y does not apply). Model VI differs from 
model C by having an exponential viscosity increase of a factor 
10 across the mantle depth range, in addition to temperature- 
dependent viscosity. Magnetodynamo diagnostics are onset 
time of inner core formation (tons), time of dynamo cessation 
due to core freeze-out (tmag), onse time of thermal strati- 
fication (t s trat, only applies if < t mag ), and the maximum 
intensity of the dipole field at the CMB (B™mb)- 
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Figure 1. Convective flow field (arrows) and potential 
temperature (colors) at a model time of t — 4.4 Gyr in 
a 2-D, annular model of mantle convection for a nominal 
Earth- like scenario (C in Table 3). Dimensionless val- 
ues of the mantle potential temperatures correspond to 
a range of 273-2673 K. 
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Figure 2. Statistically steady state surface heat flow for 
models with (time-independent) mantle internal heating 
H m =2, 4, and 6 xlO -12 W kg -1 (circle, diamond, and 
triangle, respectively). The model with H m — 4 x 10 -12 
W kg -1 (diamond) is reference scenario C (Figure 1, 
Table 3). The solid line shows the theoretical trend 
from boundary layer theory for internally heated systems 
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Figure 3. (a) Statistically steady state CMB heat flow, 
Qcmb, as a function of the CMB temperature, Tcmb- 
Results are for mantle convection scenarios C (cold sur- 
face), H (hot surface), A (surface with asymmetric tem- 
perature), SL (stagnant lid), and VI (depth-increasing 
viscosity) (see Table 3). Filled symbols indicate cases 
where bottom heating dominates over internal heat pro- 
duction in the mantle (U<0.5) (b) Heat flow curves re- 
cast in terms of dimensionless parameters Nucmb (Eq. 
10) and RacMB (Eq. 11). The grey line depicts the 
power-law relationship Nucmb oc Ra 1 ' 3 predicted from 
boundary layer theory. 
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Figure 4. Calculated inner core radius after a model 
time of t —4.5 Gyr as a function of initial CMB tem- 
perature (Tca/b ) an d present-day radiogenic mantle heat 
production rate {H mo ). The black contour marks Earth's 
present-day inner core radius of 1220 km. Labels indi- 
cate 3 scenarios that reproduce a 1220 km inner core at 
£=4.5 Gyr by using respective values of H mQ and Tcmb- 
(1) 3.8X10" 12 W kg" 1 , 5000 K; (2) 4.6xKT 12 W kg" 1 , 
4500 K; (3) 6.4xl0 -12 W kg" 1 , 4000 K. Corresponding 
dynamo evolutions for these cases are shown in Figure 
5d. 
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Figure 5. Thermal evolution of the core for nominal 
Earth-like scenario C. (a) Inner core radius. The purple 
circle shows Earth's present-day value, (b) Energy and 
(c) entropy terms (Eqs. (2) and (6) and Table 2). La- 
bel subscripts K, L, S, and G are for conduction, latent 
heat, secular cooling and graviational energy release, re- 
spectively. Qcmb is the CMB heat flow and $ is the 
entropy production rate associated with ohmic dissipa- 
tion. The initial CMB temperature used in this model 
is Ic/m -8=5000 K and the present-day mantle heat pro- 



duction rate is H„ 



3.8 x 10 -1 ^ W kg -1 (case 1 in 



Figure 4). The purple bar in (b) indicates a range of re- 
cent CMB heat flow estimates for Earth [Lay et al., 2008]. 
(d) Evolution of <3> for 3 models that reproduce Earth's 
present-day inner core radius for H mo and Tq^ b values 
as indicated by the corresponding points labeled (1-3) in 
Figure 4. 
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Figure 6. Dynamo evolution for 5 scenarios of man- 
tle convection (Section 2.4, Table 3). (a) Heat flow at 
the CMB, Qcmb- Time-dependent mantle heat pro- 
duct ion_is_taken into account by rescaling of the steady 
state Qcmb vs. Tcmb parameterizations shown in Fig- 
ure 3a. (b) Magnetic dipole field intensity calculated at 
the CMB, Be mb- The star and diamond indicate max- 
imum CMB field intensity, Bcmb, an d time of dynamo 
cessation, t mag , for the nominal Earth-like scenario C. 
The respective symbols correspond with those in Figure 
7. 
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Figure 7. Sensitivity of (a and c) dynamo cessation 
time tmag, and (b and d) maximum dipole field intensity 
at the CMB, Bq^ b , as a function of (a and b) initial 
CMB temperature T^j B and core radius tcmb, and (c 
and d) core heat capacity k c and solidus temperature at 
the planet's center T so i. Stars and diamonds respectively 
indicate Bq^ b and tmag for reference scenario C (Table 
3) and correspond with the symbols in Figure 6b . 



